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Abstract 

Based on boundary condition and dispersion relation, the superluminal group velocity in an 
anisotropic metamaterial (AMM) is investigated. The superluminal propagation is induced by the 
hyperbolic dispersion relation associated with the AMM. It is shown that a modulated Gaussian 
beam exhibits a superluminal group velocity which depends on the choice of incident angles and 
optical axis angles. The superluminal propagation does not violate the theory of special relativity 
because the group velocity is the velocity of the peak of the localized wave packet which does not 
carry information. It is proposed that a triglycine sulfate (TGS) crystal can be designed and the 
superluminal group velocity can be measured experimentally. 
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I. INTRODUCTION 



In 1968, Veselago firstly introduced the concept of left-handed material (LHM) in which 
both the permittivity e and the permeability fi are negative Q|. Veselago predicted that LHM 
would have unique and potentially interesting properties, such as the negative refraction 
index, the reversed Doppler shift and the backward Cerenkov radiation. LHM did not 
receive much attention as it only existed in a conceptual form. After the first experimental 



observation using a metamaterial composed of split ring resonators (SRR) 
of such materials has received increasing attention over the last few years 
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15j. As noted earlier, both e and /i are necessarily frequency dispersive 
in LHM. Since the frequency dispersion is important, the superluminal propagation in the 
LHM takes place flfl J 

While negative refraction is most easily visualized in an isotropic metamaterial, negative 
refraction can also be realized in anisotropic metamaterial (AMM), which does not neces- 

m nnn 



sarily require that all tensor elements of e and /x have negative values 
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Recently, Smith et al. have shown experimentally that an AMM slab designed to provide 
a permeability equal to —1 along the optical axis, will redirect E'-polarized electromagnetic 



waves from a nearby source to a partial focus 

In the present letter, we will investigate the superluminal group velocity of waves in an 
AMM. The superluminal propagation is induced by the hyperbolic dispersion relations as- 
sociated with the AMM. We describe a modulated Gaussian beam incident on the triglycine 
sulfate (TGS) crystal, which demonstrates in a straightforward manner that the peak of the 
localized wave packet displays interesting superluminal behavior. 

II. HYPERBOLIC DISPERSION RELATION 



For anisotropic media, one or both of the permittivity and permeability are second- 
rank tensors. We assume that the permittivity and permeability tensors are simultaneously 
diagonalizable. If we take the optical axis as the z axis, the permittivity and permeability 
tensors have the following forms: 
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FIG. 1: The circle and the hyperbola represent the dispersion relations of free space and AMM, 
respectively. The incident wave vector k is parallel to the group velocity v g in free space. Because 
of the anisotropy, v g = V q u;(q) in the AMM is not necessarily parallel to the refracted wave vector 
q. V q o;(q) must lie normal to the frequency contour, w(q) = const. 

where and are the relative permittivity and permeability constants in the principal 
coordinate system (i = x,y,z). It should be noted that the real AMM constructed by SRR 
is highly dispersive, both in spatial sense and frequency sense ^,15]. So these relative 



values are functions of the angle frequency u>. A general study on the shape of the dispersion 
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relation as function of the sign of these parameters has already been offered in Ref 
this work, we are interested in the case of AMM with hyperbolic dispersion relation. 

Without loss of generality, we assume that the wave vector locate at the x — z plane 
(k y — q y — 0). Maxwell's equations yield a scalar wave equation for //-polarized field . In 
free space, the accompanying dispersion relation has the familiar form: 

where k x and k z are the x and z components of the incident wave vector, uj is the frequency, 
and c is the speed of light in vacuum. We assume that there is an angle ip between the optical 
axis and the z axis. For the given polarization, the waves equation yield the dispersion 
relation in AMM as 
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a( lx + PQz + 7Qxq z = -w, (3) 



c 



where q x and q z represent the x and z components of refracted wave vector, a, (3 and 7 are 



given by 

a = (e x cos 2 if + e z sin 2 </?), 

(3 = (e x sin 2 ip + e z cos 2 ip) , 

7 = (e 2 sin 2ip — e x sin 2y?). (4) 

^■x^zf^y 

We show the dispersion geometry in Fig. ^ where a plane electromagnetic wave is incident 
from free space into the AMM. We choose the z axis to be normal to the interface, the x 
axis in the plane of the figure, and the y axis out of the plane of the figure. 

We assume here that the electric field is polarized along the y axis. The z-component of 
the wave vector can be found by the solution of Eq. Q, which yields 

1 



qz 2(3 



(5) 



the choice of sign of q z ensures that light power propagates away from the surface to the +z 
direction. Due to the anisotropy, the transmitted wave components may refract at slightly 
different angles. The values of refracted wave vector can be found by using the boundary 
condition and hyperbolic dispersion relation. 

We now determine the angle of phase refraction. The incident angle of light in free space 
is 9 1 = ta.n~ 1 (k x /k z ) and the refraction angle of the transmitted wave vector in AMM can 
be found by 6? = tan~ l (q x /q z ). From the boundary condition at the interface z = 0, the 
tangential components of the wave vectors must be continuous, i.e., q x = k x . Thus the 
refracted angle of wave vector in the AMM can be easily obtained. 

It is well known that the group velocity in anisotropic media differs from the direction of 
its wave vector. The group velocity is a very important physical quantity because it identifies 
the speed of the maximum intensity of wave packet. The group velocity in anisotropic media 



can be defined as 
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v g = V<Mq) = — e x + —e z , (6) 

where V q denotes the gradient of ^(q) in the wave vector space, e x and e z are unit cartesian 
vectors. Because of the anisotropy, v g = V q co , (q) is not necessarily parallel to the wave vector 

q 

The group velocity specifies the velocity of the peak of the wave group, and is not nec- 
essarily parallel to the wave vector. When the interface is aligned an angle with the optical 



axes of the AMM, the hyperbolic dispersion relations will exhibit some interesting effects. 
As shown in Fig. ^ the magnitude of q varies as a function of its direction. When the 
refractive wave vector q is approximately parallel to the asymptotic line of hyperbola(dash- 
dot-doted line), q can be very large. The wave front travels in the AMM with the velocity of 
v p = ou/q, so the ultra-slow phase velocity can be expected in the medium with hyperbolic 
dispersion relation. The waves propagate in different directions or in the AMM, Aq may 
not be parallel to q. Therefore, the group velocity is generally not parallel to the phase 
velocity. If q propagates in some special direction, Aq — > and the superluminal group 
velocity can be deduced. It should be mentioned that the interesting properties never exist 
in the medium with circle or ellipse dispersion relation. 



III. SUPERLUMINAL GROUP PROPAGATION 

The negative refraction index imposes that the SRR presents frequency dispersion 
This, in turn, conveys that actual anisotropic SRR is highly lossy. So the superluminal 
group velocity is not easily experimentally verified in SRR metamaterial. However, an 
extremely promising material, such as the TGS crystal, satisfies the conditions e x > 0, e z < 0. 
In principle, at far infrared frequencies due to polarization dependence of certain phonon 
resonances, the material with suitably low absorption exist 22, Thus the TGS slab 
can be prepared at low temperature and the superluminal group velocity can be measured 
experimentally. 

Based on the dispersion relation of Eq. Q, the group velocity can be derived as 

c(2aq x + jq z )e x + c(^q x + 2(3q z )e 
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(7) 



2 V a( lx + !Qxq z + I3q 2 z 

If the width of the distribution q(k) is small compared to the range of k over which k varies 
significantly, we can simply evaluate q(k) = q(k ) to obtain a good approximation. In Fig.|21 
the group velocity in the AMM is plotted. It should be noted that the superluminal group 
velocity in the AMM is induced by the hyperbolic dispersion relation. The superluminal 
group velocity does not means the propagation of energy will be superluminal. Using the 
new definition of electric and magnetic energy the energy velocity can not be 

superluminal and the conservation of energy is not violated. 

Note that a wave group can be formed from planes with different frequencies u or from 
plane waves with different k vectors [3]. To obtain a better physical picture of superluminal 
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FIG. 2: The superluminal group velocity in an AMM with different optical axial angles of (p = 
(solid), 7r/6 (dashed), 7r/4 (dash-doted) and 7r/3 (dash-dot-doted). The parameters of AMM are 
choose as e x = 1, e z = — 5 and \i y = 1 (TGS crystal). 

group velocity in AMM, a modulated Gaussian beam of finite width can be constructed. 
The field intensity distribution in free space is obtained by the Fourier integral and angular 
spectrum representation j^. Following the method outlined by Lu et al. [21], let us consider 
a modulated beam is incident from free space 

/+oo 
dk±f(k±) exp[i(k + k ± ) • r - iu Q t], (8) 
■oo 

where is perpendicular to k and u = ck Q . A general incident wave vector is written as 
k = ko + kj_. we assume its Gaussian weight is 

/(A; ± ) = ^exp[-a 2 A;J, (9) 



where a is the spatial extent of the incident beam. We want the modulated Gaussian beam 
to be aligned with the incident direction defined by the vector k = k cos 9ie x + k sin6 , /e z , 
which makes an angle Qj with the surface normal. 

Due to the anisotropy, the transmitted wave components may refract at slightly different 
angles. When the beam enters the AMM, it will no longer maintain Gaussian, but becomes 
a localized wave packet. Matching the boundary conditions for each k component at z = 
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FIG. 3: The modulated Gaussian beam incident from free space to AMM. The center wave number 
is ko = y/b with different incident angle 9j = —it/ 6. The optical axis angle tp = 7r/4 and spatial 
extent of the incident wave packet is a = 1. The peak of localized wave packet travels from point 
1 to point 2 during a cycle At with superluminal group velocity. The anisotropic parameters are 
same as those used in Fig. |2J 

gives the complex field in the form 



where T(k) is the transmission coefficient of the wave packet. The transmission coefficient 
can obtain a good approximation to simply evaluate this quantity at k . Then the trans- 
mission coefficient of the wave packet is simply given by T(ko). The normal component of 
refracted wave vector q z can be expanded in a Taylor series to first order in k to obtain a 
better approximation 



The distribution of the transmission field in the AMM can be derived from Eq. (fTU|) under 
the above approximation. 

A close look at the localized wave packet shows that the superluminal propagation of the 
peak is induced by the hyperbolic dispersion relation in the AMM. Fig. EJshows a closes view 
of the field intensity distribution of the wave packet propagating from the free space into 
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the AMM. In a circle At = 2-k/ujq of the modulated Gaussian beam, the peak of localized 
wave packet travels from point 1 to point 2. In Fig. El we set the center wave vector with a 
incident angles 6j = —71/ 6. We mark on each position of the peak of wave packet at each of 
the two times. In the At = 0.094ns, the peak of wave packet moves from (1.284, 2.571) to 
(3.144, 6.147). This propagating velocity corresponds to 0.427m/ns, which is almost exactly 
the analytical group velocity of 1.424c in Eq. (|7Jl. and the superluminal group velocity is 
demonstrated theoretically. 



IV. DISCUSSION AND CONCLUSION 



It should be noted that the superluminal group velocity has completely different origin 
from those described in isotropic LHM jf], Q, Isj or ultracold gas of atoms [2^ , where the 
superluminality is caused by the frequency dispersion of the medium. In the case discussed 
here, the superluminal group velocity in the medium is induced by the hyperbolic dispersion 
relation in the AMM. As far as we know, this kind of superluminal group velocity has not 
been recognized before. It should be mentioned, however, that the shape of localized wave 
packet is distorted once the modulated Gaussian beam is incident into the AMM. Our results 
do not violate causality or special relativity, because the group velocity is the velocity of the 
peak of the wave packet, which does not carry information. It is shown that the TGS crystal 
slab can be designed and the superluminal group velocity can be measured experimentally. 
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